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Abstract 

The accuracy of calculation of spectral line shapes in one-dimensional 
approximation is studied analytically in several limiting cases for ar- 
bitrary collision kernel and numerically in the rigid spheres model. It 
is shown that the deviation of the line profile is maximal in the center 
of the line in case of large perturber mass and intermediate values of 
collision frequency. For moderate masses of buffer molecules the error 
of one-dimensional approximation is found not to exceed 5%. 

1 Introduction 

Modern high resolution spectrometry of molecular gases requires a precise 
knowledge of spectral lines shapes. The present-day experimental techniques 
have reached so high precision and accuracy, that models describing line 
shapes have to take into account such fine effects as velocity dependence of 
collisional width and shift, correlation between velocity changing and phase 
shifting collisions, finite impact time, and radiation relaxation. The profile 
of an isolated spectral line is usually obtained by solving quantum kinetic 
equation for the off-diagonal element of the density matrix describing the 
active gas [22]. The essence of any model describing line shapes is the way 



it accounts for the collisions of radiator and perturber molecules. In prin- 
ciple, in the impact approximation the term accounting for collisions can 
be obtained by averaging the corresponding transition frequencies over ve- 
locities of the buffer molecules, as described in [21 J . The kinetic equation 
is a three-dimensional integral equation. Actually, it can be reduced to a 
two-dimensional, because the problem has an axial symmetry with respect 
to the wave propagation direction. Due to computational difficulties, this 
ab initio approach is rarely used for approximating experimental line shapes. 
Instead, most works on this subject utilize various simplified models of colli- 
sion term. The simplest of these models lead to one-dimensional equations. 
They are the strong collisions [T5], the weak collisions [7] and the Keilson- 
Storer [H] models. The more complicated ones, such as the rigid spheres [JD] 
and "kangaroo" [5] models yield 3D equations. 

The work [TD], in which the first attempt was made to proceed from the 
simplest models to a more complicated one introduced a combined approach. 
The rigid spheres collision kernel, which was analyzed in this work, leads to 
a 3D integral equation. However, in order to simplify the numerical calcula- 
tions, the kinetic equation was reduced to a ID integral equation by averaging 
the kernel over the transverse components of velocity with Maxwellian weight. 
This reduction is equivalent to adopting the assumption that the distribu- 
tion of the off-diagonal element of the density matrix over the transverse (to 
the wave propagation direction) velocity components is Maxwellian. In other 
words, this approximation, known as the one-dimensional model, neglects the 
transfer of nonequilibrium created by the interaction with the light wave to 
the distribution over the transverse velocities. Later the ID approximation 
was used in [21] and [3] to analyze general aspects of influence of collisions 
on spectral line profiles. The ID approximation in the rigid spheres model 
was utilized in [TT] to fit the experimental line shapes. The line shapes in ID 
rigid spheres were comprehensively studied and compared to other models in 
[23]. 

In our time the calculation of a line shape in the rigid spheres model can 
be performed on a usual desktop computer without utilizing ID approxima- 
tion. However, it may still become useful when improvement of precision of 
experimental measuring of spectral lines profiles will persuade researchers to 
turn to even more realistic models, with kernels describing simultaneously 
dephasing and velocity- changing collisions. 

The investigation of precision of the ID approximation is also motivated 
by the following. It has been pointed out recently [21] that probably much 
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of the disagreement between theoretical and experimentally measured line 
shapes could be removed if the calculation of dephasing term was performed 
using the correct velocity distribution of the off-diagonal element of the den- 
sity matrix instead of Maxwell's distribution. The understanding of to what 
extent this distribution really differs from equilibrium is required to clarify 
this problem. The answer to the question, how much of the non-equilibrium 
is transported to the distribution over transverse velocities, could contribute 
to such understanding. 

The question of precision of the ID approximation was first studied quan- 
titatively in [19]. In this article the problem was studied in the so-called 
"kangaroo" model. The accuracy of the ID approximation was found to be 
good. Later the same result was extended to the accuracy of ID approxima- 
tion in describing the light induced drift effect [16] in the "kangaroo" model. 
This research was continued in [17] with special attention to Dicke narrowing 
effect and gave the same result. The overall conclusion of these three articles 
regarding ID approximation is, its accuracy is high and it can be used for 
studying a wide range of problems both in linear and nonlinear spectroscopy. 
This thorough analysis does not seem quite general and comprehensive. The 
mere fact that the three mentioned simple models of collision integral lead 
to ID kinetic equation (and so the ID approximation is precise for them) 
shows that the accuracy of this approximation is determined by some fine 
properties of the collision integral. Thus the model of collision integral which 
is used to study this problem should be more realistic than the degenerate 
(in the sense that it has infinite number of eigenvectors with zero eigenvalue) 
collision kernel of the "kangaroo" model. It seems more appropriate to uti- 
lize the rigid spheres model for this purpose. This model's collision kernel is 
obtained by direct averaging of the cross-section (though for a not very real- 
istic potential), so it has many realistic features. Another factor that draws 
attention to this model is that it is often used (in combination with other 
terms accounting for collisional broadening and shifting) to fit experimentally 
obtained profiles, e.g. in [25| HI]. 

The goal of the present work is to study the precision and the area of ap- 
plicability of the ID approximation quantitatively in the rigid spheres model. 
This paper is organized as follows. In the next section the rigid spheres model 
is briefly described and the ID approximation is introduced. Section [3] is de- 
voted to analysis of the accuracy of ID model in several limiting cases. The 
results of numerical calculations in the rigid spheres model are presented and 
discussed in Section HI The conclusions are drawn in the last section. 
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2 Basic equations 



2.1 The rigid spheres model 

Let us briefly remind the basic equations defining spectral line shapes, fol- 
lowing [21]. The gas interacting with radiation can be described in terms of 
density matrix p(v), depending on velocity v. Let the wave's frequency be 
close to resonance with the transition between the levels m and n. Then, in 
the resonance approximation, the line shape is defined by the off-diagonal el- 
ement Pmniy). Below we will denote it as p(v) without indices. It is governed 
by the following master equation: 

(-no + ikv + T) p(v, u) = W(v) + Sp(v, u), (1) 

here u is the detuning of the wave's frequency from the resonance, k is 
the wave vector of the light wave, VT(v) is Maxwell's distribution, V is the 
relaxation constant, S is the collision operator. The line shape is given by 
the formula 

I(cu) = - Re J p(v) d\. (2) 

The relaxation constant V is assumed to be velocity- independent. In this 
case, it influences the line shape in the following way: 

J( M ,r) = //fr/,o) F + v- < 3 > 

that is, the line shape is a convolution of the form of the line in the absence 
of r with a Lorentzian profile. Thus the relaxation tends to conceal any 
details of the collision integral, including the differences between the initial 
3D collision integral and its ID analog. In this paper we consider equation ([1]) 
only with T = 0, because in this case the inaccuracy of the ID approximation 
is most noticeable. 

We consider the collision operator S in the impact approximation and 
also assume that the cross section of scattering of active molecule by buffer 
gas molecules is independent of the molecule's state (full phase memory). 
Under these assumptions the collision integral takes the form 

Sp(v) = / A(v, v')p(v')dV - Kv)p(v). (4) 
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The function A is the collision kernel, v is the scattering-out frequency. Ow- 
ing to the assumption we made, the functions v and A are bound by the 
following relation: 

„(v) = I A{V,v) dV. (5) 



(6) 



It can be inferred from equations (jl]) and ([5]) that 

dv S p = 

for any p(v). Another general property of the collision integral 
derived from the detailed balancing principle: 

A(v,v')W(V) = A(V,v)W(v). 

Taking into account the definition §5§ of v we obtain 

SW(v) = 0. 



can be 



(7) 



(8) 



Below we utilize the rigid spheres (or "billiard balls") collision kernel 
in all the calculations requiring a definite collision integral. In this model 
the differential cross section of scattering is considered to be independent 
of relative velocity and equal to da /do = a 2 /4, here a is the effective sum 
of radii of active and buffer molecules, and do is the element of solid angle. 
The collision kernel in the rigid spheres model is obtained by averaging the 
probability of a collision, changing the velocity of active molecule from v to 
v', over the velocities of buffer molecules. Its explicit form is [2TI [TP] 



A RS (v|v') 



N b v bT a 2 



A 2 C 



exp 




(9) 



here m is the mass of the buffer molecule, M is the mass of the active 
molecule, £ = v — v', A = |fv br , v bT = v T / yf]3 is the most probable velocity 
of the buffer molecule, f3 = m/M is the mass ratio, fi = mM jim + M) is 
the reduced mass. The kernel (Q has a singularity l/( caused by the energy 
conservation. The absence of such singularity in most phenomenological 
kernels corresponds to suppression of small angle scattering. The kernel 
(Q explicitly depends on (3 and demonstrates correct behavior in the limits 
P — > and (3 — > oo. 
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The scattering-out term in this model has the form: 



N b v bT a 2 



+ 



7T 



2z + - 

z 



Erf(z) 



where Erf [z] = J e a ' 2 <ix denotes the error function [2 J and z = v /v^x- 



2.2 The one-dimensional model 

Equation flTJ is a two-dimensional (due to axial symmetry of the problem) 
integral equation. Its direct numerical solution presents certain difficulties, 
so various simplified models are widely used. One of such models is so-called 
one-dimensional approximation. In this models the dependence of p(v) on 
the transverse components of the velocity v is considered to be Maxwellian. 
This assumption allows to write instead of (00) an equation for dependance 
of p on the longitudinal velocity v z : 

(— ko + ikv z + T) p(v 2 ) = 



v 1D p(v z ) + A 1D (v z \v z )p(v' z )dv' z , (10) 



where 



and 



A(v|v 



-v?/4 



7T 



-d 2 v ± d 2 v 



v 1D (v, 



A iD (v z ,v z )dv:. 



(11) 



(12) 



Here A 1D and v 1D are one-dimensional collision kernel and scattering-ou 
t frequency, vi and are the components of velocity, orthogonal to the 
wave propagation direction, and v^ is the most probable velocity of active 
molecules vt = v2mT. 

In the ID approximation the rigid spheres collision kernel (j9]) is reduced 

to 



7cN b a y t 



4(3 



1 + aErf (- ^ 



1 ,P 
v z + 



+1 - aErf 



2v^ 

v. 



(13) 
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(a) (b) (c) 

Figure 1: Shape of the "symmetrized" kernel A 1D (v z , v' z ) exp [(v 2 . — v^ 2 ) /2] 
at different (3 = 1 /5(a), f (b), 5(c). 

where a = sign(v 2 — v' z ). This kernel as a function of the initial and final 
velocity is presented in fig. [TJ In order to make this figure more illustrative 
we plot "symmetrized" kernel A 1D (v z , v' z )exp [(v 2 — vf) /2]. This function 
is a symmetric function with respect to transformation v z <-> v' z . At small 



perturber mass ^ < 1, Fig. 1(a), the kernel function has sharp peak near 
v' z = v z , which means that the small velocity change is the most probable. 
For comparable perturber (3 ~ 1, the peak broadens. For heavy perturber 



(3 1, Fig. 1(c), the additional ridge near v' z = —v z arises corresponding to 
elastic backward scattering on a perturbing molecule. 

Integrating expression ffl3l) over v z , we obtain the out-scattering fre- 
quency of ID rigid spheres model: 



v z Erf(V^) + — + e v " J Erf(v^t) e"* 2 dt 

(14) ' 

The line shape in the ID rigid spheres model can differ significantly from its 
shape in full rigid spheres model, as it is demonstrated in Fig. [2j It can be 
seen on this picture, that the ID line almost coincides with 3D in the wings 
of the line, but goes considerably lower in its center. In the next section we 
will analyze this difference in several limiting cases. 

The essence of the ID model is the assumption that the transfer of non- 
equilibrium to the distribution of p over the transverse components of velocity 
is negligible. Figure [3] demonstrates, that this transfer is not weak in general 
case. Fig. [3] (a) represents the distribution of p(v), divided by Maxwellian 
distribution, in absence of collisions. It can be seen that the dependence on 
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Figure 2: Shape of the spectral line I(uj) in the rigid spheres model at 
P = 100, v d = 3kv T . The solid line corresponds to the rigid spheres ap- 
proximation, the dashed line corresponds to its ID approximation. 




0.5 V Z 



(a) 



(b) 



Figure 3: Distribution of the value p(v)/W(v) over longitudinal and trans- 
verse components of velocity in the rigid spheres model at ui = and 
T = 0.5 lev?- The first plot (a) corresponds to absence of collisions. The 
second plot (b) corresponds to the case (3 — > oo, C = 3 (see section 13^]) . 
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transverse velocity is uniform. Fig. [3] (b) represents the distribution of the 
same value in presence of the collision integral. The dependence on the trans- 
verse velocity becomes strongly nonuniform. Thus it could be expected that 
the impact of transfer of non-equilibrium on line shapes would be also strong. 
However, as it will be shown, in most cases it turns out to be numerically 
small. 



3 The limiting cases 

3.1 The wings of spectral line profile 

Let us at first consider the spectral wings of the line. Note that the calcula- 
tions below take into account only rigid spheres elastic collisions and ignore 
other effects determining the far wings. It is convenient to turn to the time 
domain for this purpose and to introduce the Fourier transform of p(v, to) 

P(y,t) = f p(v,u,)e-^. (15) 
The function p(v, t) satisfies the following evolution equation: 

^ + ikv) p(v, t) = Sp(v, t) + W(y)S(t). (16) 

This equation has only one solution that tends to zero at t — > ±oo. Let us 
define the autocorrelation function \l/(i) as 

*(*) = f p(v,t)dv. (17) 



This function is different from usually used function $(i) [21] in the following 
way: \& is an even function of t, while $ is zero at negative values of t. It is 
evident from ffT5]) . ffTBT) and (fT7|) that ^(t) is a real continuous even function 
connected with the line shape by the Fourier transform 

I{u) = I J m(t)e iu)t dt. (18) 

The asymptote of I(u) at uj — >• oo is determined by the discontinuities of 
the derivatives of ^{t) at t — 0. Using (fT5"l) and taking into account (jS]) and 
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(JHD we find that (in the absence of T) the lowest order of derivative of *&(t) 
having a jump at t = is third. Thus the asymptote of I{uj) can be written 

as 

I(u) = — [ dvkv S kv W + 0(l/u 6 ) atcj^ioo. (19) 

7TCJ 4 J 

In order to obtain the value of the coefficient of the dominant term of this 
asymptote we have to substitute the definition (HI) of S into Eq. (fl9l) with some 
specific kernel A(v, v'). In the ID approximation we would have to substitute 
the reduced kernel ffTTl) . It is clear that the expression in ID case is the same 
as in initial model, the only difference is the order of integration. Thus 
the ID approximation always gives correct principal term of the asymptotic 
expansion of the line shape at u — > oo. It should be noted that if V is nonzero, 
the tails of the line have universal Lorentzian form V / (ttuj 2 ) regardless of any 
details of the collision integral. 

The fact that the tails of the line profile are insensitive to the transfer of 
non-equilibrium to the distribution over the transverse velocities means that 
the deviation of the line shape in ID approximation manifests itself mostly 
close to the center of the line. This tendency is evident, for instance, in 
Fig. El where the difference between ID and 3D line shapes is largest in the 
center. So in the rest of this paper we will mainly analyze the behavior of 
intensity in the center of the line I a = I(u = 0). 



3.2 Small collision frequency 



The next limiting case we would like to consider is the case of low collision 
frequency. In absence of the collision integral and relaxation the line shape 
is Gaussian near the center: 



exp 



(kv T y 



+ 0(u). 



(20) 



If the collision frequency is small, the line shape I(uS) can be expanded 
in power series with respect to this parameter. The first term of this series 
can be symbolically written as 



dv- 



1UJ 



lkv mjj — lkv 



W. 



(21) 



Since S is a linear integral operator, the expression ( [2T]) contains inte- 
grations over both initial and final velocities. If we perform the integration 
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only over the transverse velocities, we obtain the first term of this expansion 
corresponding to ID collision integral. Hence it is evident that the first terms 
of the power expansion in S in ID approximation and in initial model always 
coincide. 



3.3 Small perturber mass j3 <C 1 

In this case the weak collision model is applicable. In this model the colli- 
sion integral is replaced by a Fokker-Planck differential operator, and master 
equation takes the form: 

(-no + ikv) p = v ^A + Vvj p. (22) 

It is obvious that in this equation variables can be separated, and thus the 
assumption of ID approximation holds true. The collision operator of the 
weak collisions model is really the leading term of expansion of collision 
integral in power series in small mass ratio j3. In this section we are going to 
analyze validity of ID approximation in the next-to-leading order in (3. 

This analysis should exactly account for collision frequency z/, because it 
is not supposed to be small, although it is proportional to (3. For instance, 
for the rigid spheres model 

u=^N b v bT a 2 p(l+P) 2 . (23) 

Generally, the weak collisions model is applicable when (5 is small, but N b is 
large, so that v is finite. 

The term of the next-to-leading order in the collision operator is a fourth 
order differential operator. For the rigid spheres it has the form 



L = (3v 



ir ( 3 + «V+ ( 2 + j^u 2 } A+ 



+\u,iU k ViV k + ^hVA + ^A 2 

O 



(24) 



where u = v/v T and all the differentiations are with respect to u. In this ex- 
pression, summation over repeating indices is assumed. The next-to-leading 
term in the spectral line shape (that is, exact in v and of the first order 
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in (3) can be calculated by considering this operator as a perturbation and 
utilizing the time domain Green's function of equation (1221) . But here we are 
interested only in the difference of such terms in an arbitrary model and its 
ID analog. This difference can be symbolically written as 

t 

5¥ l) {t) = j dt' J dw G{t - t') (L - L m ) G{t')W. (25) 
o 

Here G(t) is the time domain Green's function of (122]) . its velocity arguments 
are omitted, and L 1D is one-dimensional fourth order differential operator 
representing one-dimensional collision operator in the considered approxi- 
mation. Let us consider the expression (|2"5]) consequently. The expression 
G(t')W (depending on velocity v and time t') is a product of Maxwellian 
distribution over the transverse velocities and a non-equilibrium distribution 
over the longitudinal velocity v z . According to the definition (jTTI) . the inte- 
gral over the transverse velocity f d 2 v±L 1D G(t')W is equal to its 3D analog 
for any value of v z : 

J d\ ± G(t - t') (L - L m ) G{t') W = 0. (26) 

The Green's function G conserves the property of having zero integral over 
for all v z , so the expression (|2"5"|) vanishes for all t. It means that ID 
model gives correct result for the correction of the first order in (3 to the 
weak collisions model. 



3.4 Heavy perturber (3 ^> 1 

The last limiting case we are going to consider is the Lorentz limit (3 — > oo 
(see [20J) of the rigid spheres model. In this case the collision kernel and 
frequency are 

A(v,v') -> ^N b a 2 5(v 2 - v' 2 ), i/(v) = 7riV b a 2 |v|- (27) 

In this case we can solve equation (CE]) for uj = and find the intensity in the 
center of the line: 

2 arcctg(C) = 7rN b a 2 

° 7T 3 / 2 kv T l-Carcctg(C) ' k 1 } 
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Constant C, being proportional to number density N b of the buffer particles, 
differs from the collision frequency only by a multiplier. Since our aim is 
only to compare 3D and ID models, we do not specify this factor in current 
section. 

At large values of C Dicke effect takes place, and the asymptotic form of 
1(0) is linear in collision frequency: 

tiIo = -4^— + Oil/C) ~ 3.38-^-. (29) 

The ID collision kernel in the Lorentz limit is 

A(v z \v' z ) = ^N b a 2 [l-B (v z - v' z ) (l - e-(^ 2 -< 2 ))] . (30) 

The out-scattering frequency in this approximation is 



v(v z ) = nN b a 2 \v z + ^-e v * (1 - Eri(v z )) ) . (31) 

The problem of finding I in this case can be reduced to solving a second 
order ODE 

* e *(*£ + v )f = -2Cv.e4fM. (32) 
dv z \ v J dv z 

with boundary conditions /(oo) = and /(0) = 1. Then the intensity in 
resonance is given by 

WI ^_fr^M dv Y\ (33) 

\Jo v{y x ) J 
At large C this expression takes the form 

nl = ( [ W(v z )-^- dv z ) C + 0(1/C) ~ 2.84-^-. (34) 

The dependance of I Q on C obtained by numerical solution of fl32|) to- 
gether with dependance (1281) is presented on Fig. HI Both functions increase 
monotonously with collision frequency, approaching their linear asymptotes 
( l29l and ( |34l) . Comparing ( f29l) and ( 1341) we find that in the Lorentz limit 
the relative deviation of I D in the ID approximation is approximately 0.2. 
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Figure 4: The dependance of intensity I Q in the center of the line on collision 
frequency C in the Lorentz limit. The upper line corresponds to 3D Lorenz 
case, the lower — to the ID approximation. The dashed lines represent linear 
asymptotes (J3U) and (J2HD 

4 Numerical results 

The difference between the rigid spheres model and ID approximation is best 
characterized by dependence of I Q on the collision frequency. Before calcu- 
lating this dependence we have to choose a value characterizing the collision 
frequency. In the Lorentz limit we used C to characterize it. In general case, 
it is customary to use the diffusion frequency, defined as Vd = vt 2 /2-D , where 
D is the mutual diffusion coefficient. The value of D depends on the collision 
kernel. Which kernel should we use to define D and z/^? It is clear that the ID 
approximation should be compared with the rigid spheres at the same values 
of physical parameters, such as and a. Using its own diffusion coefficient 
for each model would break this requirement. On the other hand, using the 
rigid spheres diffusion coefficient for both models does not seem consequent. 
Taking all that into account we decided to characterize collision frequency 
by I'd — v t/2D', where D' is the first order Chapman- Enskog diffusion co- 
efficient. The calculation of D' only accounts for longitudinal motion |13j . 
Thus the ID approximation leads to correct value of D' for any kernel. For 
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the rigid spheres we have 



3 Vr h + p 

In order to solve equations ([I]) with collision kernel (Q we used a method 
based on decomposition of p(v, u) into a linear combination of Burnett func- 
tions [T2"j E]. The problem is axially symmetric, and can be reduced to two- 
dimensional by turning to variables x = v/v a T, y = cos(k, v). The Burnett 
functions have the form: 

<Ux, y) = N, nl x l L l ^' 1 / 2 {x 2 )Pi{y), (36) 
at _ / ^n!(2f + l) 

Nnl ~ V 2r(n + / + s/2)' (37) 

^ o m!(n-m)!r(m + / + 3/2 ' 1 j 

, N 1 ^ (-l) fc (2/-2A;)! , 2 , 

Here A^ n ; is normalizing factor, L^ +1//2 (x 2 ) are generalized Laguerre polyno- 
mials, P\{y) are Legendre polynomials, T(. . . ) is Euler's Gamma function [2]. 
The desired absorbtion intensity is given by n = 0, I = coefficient in the 
decomposition of p(v,u). 

For numerical solution, we need to limit the decomposition by some 
I = 'max and n = n max . The structure of emergent system of linear alge- 
braic equations allows to perform the calculation considerably faster than 
for a generic linear system. The matrix of the system turns out to be block 
tridiagonal if we group the coefficients of decomposition of p(v, u) so that I 
numerates blocks, and n numerates the elements within each block. Then the 
blocks in the main diagonal are square symmetric matrices, and the blocks in 
the neighboring diagonals are two-diagonal matrices, the blocks above and 
below the main diagonal are transposed with respect to each other. The 
vector in the right hand side of the system has only one non-zero element, 
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R = 1. Thus the system has the form 









R = 





A = 




VJ 





p = 

( Ao 


V : 



R, 

D 
Ai 





D x 
A 2 






D 2 



(40) 



(41) 



This system was solved by the special method for block tridiagonal ma- 
trix [8j. As in the scalar marching method for a tridiagonal matrix |18j . two 
sequences of coefficients are calculated using recurrent formulas. The differ- 
ence is, these coefficients are not numbers but matrices and vectors. In our 
case, only matrix sequence of coefficients needs to be calculated because the 
vector coefficients are defined by the right hand side and thus they are equal 
to zero: 

Mi — - (Ai + DiMi+iDjy 1 , 



M, 



-A 



-i 



So, initial system is reduced to a system of dimensionality n Ti 

(A + DqMiDq) po = Ro- 



+ 1: 



This system was solved using the Givens rotation method [8J. 

The one-dimensional problem fllOp was solved numerically in a straight- 
forward way by discretization of velocity. The integrals appearing in ffTOl) 
were approximated with sums by Simpson's formula with accuracy 0(1/N A ) 
[j~8] . where N is number of points. The obtained system of linear equations 
was solved using Gaussian method. 

The results of 3D and ID numerical calculations are presented on Fig. 
O [3 The lower panels show the difference between ID and 3D calculations. 
The deviations are small, the corresponding scaling coefficients are indicated. 
Fig. [5] corresponds to (3 = 1/5, fig. [B]to (3 = 1 and fig. [7] to (3 = 5. On the 
first two figures the three plots were calculated at u d /k\ T = 1/3, 1, 3, from 
left to right. On the last figure the left plot corresponds to z/ d //cv T = 1, and 
the right one to v^jkvT = 3. All profiles were calculated at 7 = 0.03. The 
3D plots were computed with l max = n max = 48. 

The typical magnitude of deviation of ID plots on the three plots of fig. 
[5]is 10" 5 , 10" 4 , 10" 3 , from left to right; on fig. E it is 10" 4 , 10" 3 , 10" 2 ; and 
on fig. [7] it is 10~ 2 and 10 _1 . It can be seen that the difference between line 
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Figure 5: The line profiles in 3D and ID rigid spheres model at (3 = 1/5 and 
v d = l/3(a), 1(6), 3(c). 
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Figure 7: The line profiles in 3D and ID rigid spheres model at (5 = 5 and 
i/ d = l(o),3(6). 
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profiles in 3D and ID case increases with (3 and with v&. In agreement with 
theoretical expectations, the deviation tends to zero at the wings of the line 
profiles. 

Taking into account the behavior of deviation of ID approximation in the 
limiting cases, it seems natural to expect that the relative deviation would 
grow monotonously with collision frequency and with (3, reaching maximum 
values of ~ 0.2. However, this assumption proves to be wrong. Numerical 
calculations in the rigid spheres model show that at (3 < 100 the relative 
deviation of ID approximation is less than 0.1. At (3 < 30 the deviation is 
less than 0.05. 

Fig. [S] presents I {yd) obtained by numerical calculation in the rigid 
spheres model and in its ID approximation for (3 = 100 in comparison with 
the Lorentz case. All the curves on this plot increase monotonously and 
approach corresponding linear asymptotes as collision frequency tends to in- 
finity. The 3D j3 — 100 curve is close to 3D Lorentz line at small v&. As the 
collision frequency increases, it goes down and crosses the ID line, represent- 
ing both Lorentz and (3 = 100 case. 

The behavior of the curves on this plot can be described as follows. The 
Lorentz model is the limiting case of the rigid spheres model at (3 — > oo. The 
dependence of I D on Vd in the rigid spheres model approaches that dependence 
in the Lorentz limit as (3 goes to infinity. However, this approaching is not 
uniform in v^'- at any fixed f3, no matter how large, at sufficiently large v& the 
dependence deflects from Lorentzian asymptote downwards and approaches 
some other linear asymptote. This other (final) asymptote goes lower than 
( I34p . so the 3D line crosses the one corresponding to ID approximation. It 
can be said that due to this intersection of the two plots the error of ID 
approximation at reasonable values of (3 is much less than in Lorentz limit. 

As it is known, in case of full phase memory in the limit of large collision 
frequency the line shape takes the form 

1 Dk 2 

'<"> = n «^ + °< D >- < 42 > 

In case of large (3 this formula gives 

rf -S(iV° (1/ ^ 177 (^ (43) 

This formula describes the asymptotic behavior of the rigid spheres model 
at large but finite (3. It perfectly agrees with numerical results. However, 
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Figure 8: The dependance of intensity in the center of the line I Q on collision 
frequency v& in the 3D and ID rigid spheres model for (3 = 100 and in 
the Lorentz limit. The solid line corresponds to 3D rigid spheres model 
at (3 — 100, the dot-dash one — to its ID approximation. The dashed line 
corresponds to 3D Lorentz limit. The ID Lorentz curve is not plotted because 
it goes very close to ID with (3 = 100. 
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formula (I43p does not describe the Lorentz limit itself, its asymptote is given 
by (12T)|) . The derivation of (jUJ) implies that the only distribution that is 
turned to zero by the the collision operator is equilibrium distribution. The 
Lorentzian collision operator is degenerate in this sense, because it turns 
to zero all distributions depending only on absolute value of velocity. This 
consideration shows that in case of large but finite /3 function I (vd) should 
behave as follows: while Vd is not too large, it is close to (|2"S1) . At large values 
of Vd the function Id^d) must approach the asymptote (|43|) . This is exactly 
what numerical calculations show. The remaining question is, at what values 
of Vd does the function I {i>d) switch from (I28I) to (143]) . 

Let us consider the collision operator of the rigid spheres model. As it 
is shown in pQ, all its eigenvalues are positive, except the one corresponding 
to equilibrium distribution, which is zero. It is clear that in case of large (3 
the minimal positive eigenvalue has the order of magnitude of (3~ l . It means 
that the coefficient in 0{l/vd) term in (H21) is of the order (3. Thus the first 
term in (j42H is much greater than the second when i^/(fcvp) 2 3> (3. So the 
transition happens (and the ID plot intersects with rigid spheres plot) at 
v d ~ y/]3 kvT- This result is in good agreement with numerical calculations. 

As it was mentioned, the most significant difference between 3D and 
ID line profile is observed in case which is effectively Lorentzian: 1 <C 
v dl (^ v r) *C \Tf3~- The reason of such big deviation is that at large f3 the 
3D rigid spheres model and ID approximation behave differently: in 3D 
model, the kinetic energy of an active molecule almost does not change in 
collisions. In contrast to that, in ID approximation there is no such conser- 
vation. This "energy persistence" property of the 3D rigid spheres collision 
kernel can be interpreted in terms of the generalization of the Keilson-Storer 
model recently proposed in j3]. This model introduces two velocity persis- 
tence parameters, 7 m and j Q . They are responsible for the persistence of the 
modulus of velocity and its orientation, correspondingly. The separation of 
these two parameters is incident in the rigid spheres model. Indeed, in case of 
heavy perturber gas, the typical change of the speed of an active molecule in 
a single collision is of the order vr/3 -1 , while its orientation changes totally. 
Thus, in case of large (3 the modulus persistence parameter 7 m is close to 
unity, and the difference has the order of magnitude 

l- lm ^(3-\ (44) 

This simple estimate agrees remarkably well with the figures obtained in 
[I] by simulation: for H 2 in nitrogen ((3 ~ 14) and in argon (/3 ~ 20) the 
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estimate ([ill) gives 7 m = 0.93 and 7 m = 0.95, while the figures presented in 
[4] are 0.92 and 0.96 correspondingly. 

5 Conclusion 

We have analyzed several limiting cases for arbitrary collision kernel and 
found the following. The term 1/u 4 in the asymptote of the tails of the line 
shape, which is principal in absence of dephasing, is given correctly by the 
ID approximation. Thus the transfer of disequilibrium distribution to the 
transverse components of velocity manifests itself mostly close to the center 
of the line, and the error of ID approximation can be characterized by the 
deviation of absorption in the center of the line. We demonstrated that in 
case of small collision frequency the first correction in this parameter is also 
given correctly by ID approximation. 

In the limit of small perturber to radiator mass ratio f3 it is known that 
the line shape in ID approximation coincides with that in the initial 3D 
model. We have shown that the line shape given by ID approximation re- 
mains correct also in the next order in (3. In the opposite limiting case of 
large f3 any realistic 3D collision integral conserves the kinetic energy of ac- 
tive molecules. The ID collision integral cannot reproduce such property. We 
have shown that due to this difference, the most significant deviation of ID 
line shapes is observed in case of heavy perturber molecules. The deviation 
vanishes at small collision frequency, and reaches its maximal value in the 
hydrodynamical limit v 3> kvT- We found that in the rigid spheres model for 
infinite (3 the relative error of ID approximation reaches the value of ~ 20%. 

The case of intermediate mass ratio was analyzed numerically in the rigid 
spheres model. The inaccuracy of ID approximation was found considerably 
smaller than expected value of ~ 20%. This happens due to intersection of 
the plots describing the collision frequency dependance of intensity in the 
center of the line in ID and 3D models. For moderate values of j3, that is, 
P Pa ~ 30, the relative error of ID approximation does not exceed 0.05. 
The error increases monotonously with (3. At P < Pq the error takes its 
maximal value in the hydrodynamical limit. At larger (3 the maximal error 
occurs at v ~ \[pkv T . 

Thus the one-dimension model could be applied for light perturbers, at 
low pressure or in the problems where 5% is sufficient accuracy. For arbi- 
trary buffer particles or precision calculation the three dimension billiard ball 
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approximation becomes preferable. 
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